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We present a phenomenological model for the remanent magnetization at low tempera- 
tures in the quasi-one-dimensional dilute antiferromagnets (CHsNH^Mni-^CdxCla^H^O and 
(CHs^NEbMni-xCdrCls^EhO. The model assumes the existence of uncompensated magnetic mo- 
ments induced in the odd-sized segments generated along the Mn 2+ chains upon dilution. These 
moments are further assumed to correlate ferromagnetically after removal of a cooling field. Using 
a (mean-field) linear-chain approximation and reasonable set of model parameters, we are able to 
reproduce the approximate linear temperature dependence observed for the remanent magnetization 
in the real compounds. 
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I. INTRODUCTION 

At low temperatures, quasi-one-dimensional mag- 
netic systems exhibit a wealth of-interesting behav- 
ior, such as dimensional crossover paa random quan-i 
turn paramagnetismJjLorder-by-disorder phenomenapEl 
and Griffiths phases pB which have motivated many ex- 
perimental and theoretical investigations. In most of 
these systems, three-dimensional (3D) order is eventu- 
ally induced by interchain interactions. Taking advan- 
tage of the many analytical results available for one- 
dimensional model systems, this situation has been mod- 
eled in a variety of ways. Most of the existen±|-ap=. 
proaches are based on linear-chain approximations Jjtt3 Ell 
which treat correlations along the chains in an exact way, 
while introducing interchain couplings via effective fields. 
These approximations have been quite successfully ap- 
plied to pure systems, andJiswe given rise to generalized 
Ginzburg-Landau theories Jala which account for fluctua- 
tions. AlsOp-th^s-Jj^y^ikeen widely used to explain disor- 
der effects, EJCJlla'OtZl which are among the main topics 
of research on quasi-one-dimensional systems. 

In the present work we consider a class of quasi- 
one-dimensional compounds,EM2l represented by 
(CH 3 NH 3 )MnCl 3 -2H 2 (abbreviated MMC) and 
(CH 3 ) 2 NH 2 MnCl 3 -2H 2 (abbreviated DMMC), consist- 
ing of localized spin systems in which the Mn 2+ ions (spin 
S = §) lie along the crystalline b axis, forming chains, 
and are antiferromagnetically coupled to each other by 
an intrachain interaction J/ks around 3 KpJvfagnetic 
susceptibility and specific heat measurements! 2 ^ indicate 
the onset of 3D long-range order at Neel temperatures 
T N = 4.12 K for MMC and T N = 6.36 K for DMMC, 
with the magnetic moments aligning along the a axis. 
These temperatures are compatible with an interchain 
interaction \J±\ ~ |J| x 10~ 2 . The character of this 
interaction is not reported in the literature. However, 
owing to the behavior of the materials upon dilution 
with non-magnetic Cd 2+ ions, it has been suggested 
(see below) that ferromagnetic interchain couplings are 



present. At temperatures above T ~ 10 K, susceptibility 
results are well described by a one-dimensional Heisen- 
berg S = | model, but at lower temperatures anisotromz. 
effects (probably of dipolar origin) become relevant .c2l 
Calculations based on a classical anisotropic Heisenberg 
model, with parameters derived from experiments, on 
DMMC, reinforce the importance of anisotropy.Efj In 
particular, the behavior of the correlation length along 
the chains is predicted to cross over from Heisenberg-like 
to Ising-like as the temperature is lowered. 

Substitution of small amounts of non-magnetic Cd 2+ 
for Mn 2+ ions-iaduces the appearance of a remanent 
magnetizationtSt£l below Tn, when samples are cooled in 
the presence of fields of a few Oe directed along the easy 
axis. This remanent magnetization is observed to vary 
linearly with temperature, except immediately below T«_ 
where demagnetization effects seem to be relevant .113 
Moreover, there is an excess parallel susceptibility, which 
is in general associated with the existence of uncompen- 
sated magnetic moments in odeUsized segments formed 
along the chains upon dilution.Eil Apparently, the linear 
temperature dependence of the remanent magnetization 
is of universal character, as inferred from measurements!!^ 
performed on DMMC doped witk.Cd 2+ (non-magnetic) 
and Cu 2+ (S — |). ExperimentsEj performed on similar 
compounds, CsMnCl 3 -2H 2 (CMC) and CsMnBr 3 -2H 2 
(CMB) doped with Cu 2+ (for which the signs of the inter- 
chain interactions are well known) , revealed that a rema- 
nent magnetization appears in CMB, in which interchain 
couplings are ferromagnetic along one of the transversal 
directions and antiferromagnetic along the other; in con- 
trast, no net magnetization is observable in CMC, where 
all interactions are antiferromagnetic. These experimen- 
tal results, combined with the observation that some ef- 
fective ferromagnetic coupling is expected in order to sus- 
tain a net remanent magnetization, have led to the idea 
that ferromagnetic interchain interactions should also be 
present in DMMC and MMCE^ However, in the lack of 
experimental data, up to now no conclusive evidence on 
this point seems to be available. 



2 



In this paper we introduce and discuss a phenomeno- 
logical model for the low-temperature magnetic behavior 
of those compounds. By virtue of the previously dis- 
cussed anisotropy effects, we believe the qualitative as- 
pects to be captured by a S — I Ising model, which in 
the pure limit (and in the simplest case) is described by 
the Hamiltonian 



n = j 



(1) 



where J > 0, r is a lattice vector, b is the primitive 
vector along the crystalline b direction, S is a vector con- 
necting a site to its nearest neighbors in the ac plane, 
Jg = J± > if d is parallel to the a axis, and Jg = —J± 
if S is parallel to the c axis. Our approach is based on a 
linear-chain approximation, which treats the intrachain 
couplings (J) exactly, while introducing the weak inter- 
chain interactions ( J± <C J) via Curie- Weiss terms con- 
necting all spins (in such a way that a staggered effective 
field results, combining both ferro- and antiferromagnetic 
interchain interactions in a cooperative manner) . At very 
low temperatures the chains are antiferromagnetically 
ordered, with a characteristic two-sublattice structure. 
Upon dilution, a very long chain breaks into finite seg- 
ments, and uncompensated magnetic moments appear 
at the ends of odd-sized segments. On phenomenological 
grounds, we assume these moments to correlate ferro- 
magnetically, with their direction determined in the ex- 
periments by the cooling field. For each segment of spins, 
the partition function can be exactly calculated; the total 
free energy of the chain is obtained by summing the free 
energies of segments of all sizes, with prop er w eighting 
factors. Thi s procedure is detailed in Sec. IIA. Then, 
in Sec. [II Bj , we include the Curie- Weiss terms and dis- 
cuss the results of the approximation. We show that this 
approach reproduces the linear temperature dependence 
of the remanent magnetization and the existence of an 
excess susceptibility. The final section is devoted to a 
discussion and conclusions. 



describe long-range antiferromagnetic order in the pres- 
ence of interchain interactions. According to the phe- 
nomenological hypothesis that there are uncompensated 
magnetic moments pointing in a preferred direction, de- 
termined by the cooling field, we assume that the spins 
at the ends of odd-sized segments are always under the 
action of a field hi. When the field is removed, the mo- 
ments would remain uncompensated due to pinning by 
the non-magnetic impurities. For even-sized segments, 
the particular choice of a field hi at j — 1 is of no con- 
sequence, since in this case the partition functions are 
symmetric under inversion. 

As we are considering finite values of n, we must treat 
separately the cases of odd- and even-sized segments. Us- 
ing the transfer-matrix technique, we can write the par- 
tition functions as 



odd 



and 



T— Ti 



Vl 



vi 



Vl ) = ( V 2 



T 2 T" 



Vl 



(3) 



(4) 



where n is an even number, T = T1T2, the elements of 
the 6x6 matrices Ti and T 2 are given by 

Ti (5,, S-) = e -l 3 JS,s ] +^p(h 1 s l +h 2 s 3 )+^pD(s 2 % +s 2 J ) ^ ^ 
T 2 (S i ,S j )=Ti(S j ,S i ) ) (6) 



and the components of the vectors Vi and v 2 are 
Vl ( Sj ) = eiK^+DSj)^ 
v 2 ( Sj ) = e MWS3). 



(7) 
(«) 



The free energies associated with odd- and even-sized 
segments are 



Fl 



odd 



-fc B TlnZ°li 



(9) 



and 



II. PHENOMENOLOGICAL MODEL 
A. Nearest-neighbor interactions 

We initially consider an open segment of n Ising spins 
with antiferromagnetic couplings and alternating fields, 
described by the Hamiltonian 

n—1 n n 

H n = JY, SjS j+ i - Yl h & ~ D T, s l ( 2 ) 

3=1 3=1 3=1 

where J > and hj — hi (h 2 ) for odd (even) j; a crystal 
field D is also introduced. The spin variables Sj take the 
values ±5, ±| and ±|. The alternating fields are intro- 
duced to give room to a staggered effective field needed to 



F° von = -k B T\nZ™ cn . (10) 

We now consider a very long chain, and assume that 
each of the N sites is occupied by a spin with probability 
p. For < p < 1, the chain is composed of finite segments 
of spins separated by empty sites. In the — » 00 limit, 
the number of segments of size n is NP(n) — N(l — 
p) 2 p n . Assuming that each segment is described by the 
Hamiltonian in Eq. (||), the total free energy per spin is 
given by the infinite series 

fnn(hl,h 2 ,T) = IJ2 t P ( n ~ ^n-l + ^H^™ 1 ] • 
n even 

(ii) 

For p < 1, as nP(n) becomes negligible for sufficiently 
large n, this infinite series can be truncated and readily 
evaluated numerically. 
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Let us denote by type 1 (type 2) those spins under 
the action of a field hi (h-z). The numbers N% and N 2 
of spins of either type in a chain can be determined by 
noting that in a segment of size n there are n/2 type-1 
spins if n is even and (n + l)/2 type-1 spins if n is odd. 
Thus, the fractions of type-1 and type-2 spins are 



N 



E p ^ 



odd 



and 



N2 
N 



E p ^ 



n- 1 



E p^)- 2 

n even 



E^( 



n p 



odd 



l+p' 



(12) 



(13) 



respectively. For p < 1, the difference between the two 
fractions will obviously generate a non-zero magnetiza- 
tion at zero temperature. 



B. Linear-chain approximation 

In order to mimic the weak interchain coupling in the 
real compounds, we now assume that, in addition to the 
nearest-neighbor couplings inside each segment, there are 
also ferromagnetic Curie- Weiss (CW) interactions con- 
necting all spins in the chain. We further assume that 
the CW interactions between two type-1 or two type-2 
spins have strength J cw /N, but that the CW interac- 
tions between a type-1 and a type-2 spin are weaker by a 
factor 7. This 7 factor is introduced to allow for off-plane 
interchain couplings; in the pure limit (p — 1) the chains 
are expected to exhibit antiferromagnetic order, so that 
7 must be smaller than unity. Upon dilution, we expect 
the antiferromagnetic arrangement to survive inside each 
segment, and in principle this could lead to a variation of 
7 with the concentration p, since the magnetic arrange- 
ment in the planes perpendicular to the chains could be 
disturbed. In any case, our results suggest that 7 is very 
small, if not zero, in the compounds under consideration. 

The contribution of the CW interactions to the total 
energy per spin is 



-pJ C vj{rn\ + 2jmim,2 + ml), 



(14) 



where toi (7772) is 
to spins of type 1 
averages toi and 
tion of the spins, 
variables. So, we 
per spin a nn (TOi, 
interactions, and 



the magnetization per magnetic ion due 
(type 2). Since e C w depends only on the 
?Ti2, and not on the detailed configura- 
it is convenient to perform a change of 
now introduce the Helmholtz potential 
77i2, T), related to the nearest-neighbor 
defined by the Legendre transform 



a nn (mi,m 2 ,T) = f nn (hi,h 2 ,T) + mjli + m 2 h 2 , (15) 



where hi and h 2 are effective fields and 

fdfn_ 



mi 



V dhi J~ h2 , 



and 7772 



dfnn 

dh 2 / hi,T 
(16) 



For given values of mi and 7772 we can write the total 
Helmholtz potential, 



a(mi,777 2 ,T) = a nn (mi,7772,T) + e c 



(17) 



from which we obtain the relation between the external 
magnetic fields hi, h 2 and the effective fields, 



hi 



da 



and similarly 

da 

dm-. 



= hi - 2pJ cw (mi +77772), (18) 



= h 2 ~ 2pJ cw ('jm 1 + m 2 ). (19) 



2 / mi,T 

Comparing these last results (for 7 = 0) with the local 
field at a site r due to its q± nearest-neighbors in adja- 
cent chains, as given by the Hamiltonian in Eq. (Q), we 
conclude that J cw can be estimated as 



Jew — -pq±J±- 



(20) 



The thermodynamically stable magnetizations are 
those which minimize the free-energy functional 

$ (hi, h 2 , T; mi, m 2 ) = a(mi,m 2 ,T) — mihi - m 2 h 2 
= /„„(%, /72,T) -£cw (21) 

For low temperatures and small ratios J cw /J, setting 
hi = h 2 = 0, the stable values of mi and 7772 have oppo- 
site signs. In the presence of dilution (p < 1), since we 
have |mi| ^ \m 2 \, the model predicts a remanent mag- 
netization per lattice site, 777 r , given by 

m r — p(mi + m 2 ). (22) 

In the T — > limit, m r reaches a saturation value, 

Ni - N 2 p(l-p) a 
lim m r = — S = — r 6, (23) 

where S = |. The zero-field differential susceptibility Xo 
can be calculated by setting hi = hi = h and taking the 
h — » limit, 



Xo = lim 



C?777 r 



h^o dh 



(24) 



The Neel temperature is obtained from the solution of 
the equation 



d 2 $ «9 2 $ / <9 2 $ 



if m\ \dmidm 2 / 



= 0, 



(25) 



mi— 7712— 



in the absence of an external field. |— ■ 

In Fig. |l| we show the experimental datali^ for the tem- 
perature dependence of the remanent magnetization in 
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FIG. 1: Experimental data (circles) and theoretical calcu- 
lation (solid curve) for the temperature dependence of the 
remanent magnetization in DMMC with 4.5% Cd. The mag- 
netization is normalized to its value at the lowest temperature 
for which experimental data are available. 
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FIG. 2: Calculated zero-field susceptibility per magnetic ion 
in the pure limit (dashed curve) and for 4.5% dilution (solid 
curve), using the same parameters as in Fig. |lj The arrows 
indicate the corresponding Neel temperature, which is lower 
in the dilute case. The inset shows low-temperature behavior. 



DMMC doped with 4.5% Cd (this concentration was esti- 
mated from high-temperature fits to a Curie- Weiss law) . 
We also show results of our calculations for the remanent 
magnetization with 4.5% dilution, J cw /J — 1.5 x 10~ 2 , 
7 = 0, and D = 0. We obtained the best fit for the lin- 
ear portion of the curve by setting the theoretical Neel 
temperature to 1.14 times the experimental value (which 
amounts to fitting J). This is a reasonable procedure, 
since our calculations are of a mean-field nature, and thus 
we do not expect to obtain good quantitative results for 
the Neel temperature. Of course, the qualitative features 
of the calculations are not sensitive to small variations 
in the parameters; however, no strictly universal behav- 
ior (in the sense of data collapse) could be identified. We 
point out that setting the value of the crystal field to high 
positive values turns the system into a S — \ Ising model, 
and in this case the linear temperature dependence of m r 
could not be so well reproduced. Note that, in view of 
Eq. (pfj|), the value of J cw /J used in the fit is fully com- 
patible with the estimated experimental value of Ji/J 
mentioned in the Introduction. The calculated ratio of 
the Neel temperatures of the diluted and pupa systems is 
0.86, compared to the experimental estimateta of 0.99 for 
the real material. From Eq. (23), the saturation value 
of m r for 1% dilution is 0.497% of the sublattice magne- 
tization in the pure system- in excellent agreement with 
the experimental estimated of 0.5% for MMC with 1% 
Cd. 

In Fig. H we use the previous set of parameters to 
plot the calculated zero-field susceptibility xo both in 
the pure limit and for 4.5% dilution. The broad maxima 
in the curves reflect short-range intrachain antiferromag- 
netic correlations, while the cusps (indicated in the fig- 
ure by the arrows) correspond to the Neel temperatures 
(TV). As it is clear in the inset, for the dilute case we 
observe other features at lower temperatures. The small 



maximum close to T = is due to isolated spins, whose 
sole energy scale is given by the very weak interchain 
couplings, while the neighboring shoulder is due to small 
odd-sized segments, whose end spins are uncompensated 
(even-sized segments give negligible contributions to xo 
at such low temperatures). This is in sharp contrast to 
the pure limit, in which the susceptibility vanishes expo- 
nentially for T <T N . 

It should be mentioned that the present approach is 



a generalization of that used by SlotteEj to investigate 
the dilute S = \ Ising chain with competing short- and 
long-range interactions. However, owing to the presence 
of competition, his approach does not contemplate the 
possibility of long-range antiferromagnetic order at finite 
temperatures, even in the pure limit. 



III. CONCLUSIONS 

We introduced a phenomenological model for the re- 
manent magnetization (m r ) in a class of dilute quasi-one- 
dimensional antiferromagnets, composed of weakly inter- 
acting spin chains. The model assumes the existence of 
uncompensated spins at the ends of odd-sized segments 
formed along the chains upon dilution. These spins are 
supposed to remain ferromagnetically correlated after a 
cooling field is removed. By using a linear-chain approx- 
imation, in which interchain interactions are treated at 
a mean-field level, we were able to reproduce the linear 
temperature dependence of m r for a set of parameters 
compatible with experimental estimates. 

Our linear-chain approximation is based on the as- 
sumption that, even upon dilution, each segment feels 
an staggered effective field. Of course, this assumption 
is subject to some restrictions. Depending on the impu- 
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rity concentration 1 — p, the existence of uncompensated 
moments pointing in a preferred direction could lead to 
a complete destabilization of the magnetic ordering per- 
pendicular to the chains (this can be seen by considering 
the effect, in a particular chain, of two neighboring non- 
magnetic ions, which may invert the roles of the alter- 
nating sublattices). In this case, spins along the chains 
would feel the same interchain effective field, irrespective 
of their position. Actually, this would lead to a ferromag- 
netic transition (with a diverging susceptibility), and the 
long-range antiferromagnetic ordering would not be re- 
covered even as p — » 1. We have performed calculations 
near this limit, and have checked that the critical tem- 
perature depends linearly on 1 — p, being thus too small 
compared with experimental findings. Moreover, it is not 
possible to reproduce the linear temperature dependence 
of the remanent magnetization. We conclude that, at 
least in the low impurity concentrations used here, for 
which the occurrence of two neighboring non-magnetic 



atoms in the same chain is a rare event, our approxima- 
tion is reasonable. 

There remains the topic of identifying the precise 
mechanism responsible for the persistence of ferromag- 
netic correlations between the uncompensated spins. 
Monte Carlo simulations based on the Hamiltonian in 
Eq. (||) could be used to verify whether it is sufficient 
or necessary to have both ferro- and antiferromagnetic 
interchain interactions present in order to give rise to 
a remanent magnetization in quasi-one-dimensional sys- 
tems. 
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